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Abstract. Using domain wall fermions, we estimate -Bx(m ~ 2GeV) = 0.602(38) 
in quenched QCD which is consistent with previous calculations. We also find ratios 
of decay constants that are consistent with experiment, within our statistical errors. 
Our initial results indicate good scaling behavior and support expectations that 
0{a) errors are exponentially suppressed in low energy {E <C a,~^) observables. It 
is also shown that the axial current numerically satisfies the lattice analog of the 
usual continuum axial Ward identity and that the matrix element of the four quark 
operator needed for Bk exhibits excellent chiral behavior. 

1 Introduction 

We recently reported ^ on calculations using a new discretization for sim- 
ulations of QCD, domain wall fermions (DWF) ^, which preserve chiral 
symmetry on the lattice in the limit of an infinite extra 5th dimension. There 
it was demonstrated that DWF exhibit remarkable chiral behavior [Q even at 
relatively large lattice spacing and modest extent of the fifth dimension. 

In addition to retaining chiral symmetry, DWF are also "improved" in 
another important way. In the limit that the number of sites in the extra 
dimension, Ns, goes to infinity, the leading discretization error in the effec- 
tive four dimensional action for the light degrees of freedom goes like O(a^). 
This theoretical dependence is deduced from the fact that the only operators 
available to cancel 0{a) errors in the effective action are not chirally sym- 
metric ||. For finite Ns, 0{a) corrections are expected to be exponentially 
suppressed with the size of the extra fifth dimension. Our calculations for Bk 
show a weak dependence on a that is easily fit to an ansatz. Preliminary 
results for the ratios fTr/nip and fk/ fir indicate good scaling behavior as well. 

We use the boundary fermion variant of DWF developed by Shamir. For 
details, consult Kaplan]^ and Shamir[|j. See Ref. [|| for a discussion of the Ad 
chiral Ward identities (CWI) satisfied by DWF. Our simulation parameters 
are summarized in Table 1. 



2 Results 



We begin with the numerical investigation of the lattice PCAC relation. The 
CWI are satisfied exactly on any configuration since they are derived from the 
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Table 1. Summary of simulation parameters. M is the five dimensional Dirac 
fermion mass, and m is the coupling between layers s = and Ns — 1. 





size 


M 


m(# conf) 


5.85 


16^ X 32 X 14 


1.7 


0.075(34) 0.05(24) 


6.0 


16"* X 32 X 10 


1.7 


0.075(36) 0.05(39) 0.025(34) 


6.3 


24-" X 60 X 10 


1.5 


0.075(11) 0.05(15) 0.025 (22) 



corresponding operator identity. We checked this expUcitly in our simulations. 
In the asymptotic large time limit, we find for the usual PCAC relation 

which goes over to the continuum relation for aTOjr <^ 1 and Ns oo (see 
Ref. Q for operator definitions). The second term on the r.h.s. is anomalous 
and vanishes as Ng ^ oo. It is a measure of explicit chiral symmetry breaking 
induced by the finite 5th dimension. At Q/g^ = 6.0 and Ng = 10 we find the 
l.h.s. of Eqj^ to be 0.1578(2) and 0.1083(3) for m = 0.075 and 0.05, respec- 
tively. The anomalous contributions for these two masses are 2 x (0.00385(5) 
and 0.00408(12)), which appear to be roughly constant with m. Increasing 
Ns to 14 at m = 0.05, the anomalous contribution falls to (2x) 0.00152(8) 
while the l.h.s. is 0.1026(6), which shows that increasing Ng really does take 
us towards the chiral limit. 

Next we investigate the matrix element of the four quark operator Oll 
which defines Bfc. {^K\Oll\K') vanishes linearly with m in the chiral limit in 
excellent agreement with chiral perturbation theory (Fig. 1). 

In Fig. 2 we show the kaon B parameter at each value of versus 
a/jT which is used to set the lattice spacing. The results for Bk depend 
weakly on 6/5^, and are well fit to a pure quadratic in a. We find Bk^H = 
a^^) = 0.602(38) in the continuum limit. This value is already consistent 
with previous results |7| though it does not include the perturbative running 
of Bfc to a common energy scale. This requires a perturbative calculation 
to determine the scale dependence of Oll, which has not yet been donc|H. 
From Table 2, the energy scale at 6/5^ = 6.0 is roughly 2 GeV. 

At 6/g^ = 6.0, we have also calculated B^ using the partially conserved 
axial current A'^(x) (and the analogous vector current). This point split con- 
served current requires explicit factors of the gauge links to be gauge invari- 
ant. Alternatively a gauge non- invariant operator may be defined by omit- 
ting the links; the two definitions become equivalent in the continuum limit. 
Results for the gauge non-invariant operators agree within small statistical 
errors with those obtained with naive currents. Fig. 2(see Ref. |l[ for oper- 
ator definitions). The results for the gauge invariant operators are somewhat 
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Fig. 1. The matrix element {^P\Oll\P') vs. m. m is proportional to the quark mass 
in lattice units. A^^ = 10 = 6.0, 6.3) and 14 (6/fl^ = 5.85). 



larger: ''(^ = a-^) = 0.857(20) and 0.946(28) at m = 0.05 and 0.075, re- 
spectively. A similar situation holds in the Kogut-Susskind case where it was 
shown that the gauge invariant operators receive appreciable perturbative 
corrections which bring the two results into agreement |^ . 



Table 2. Lattice spacing and decay constant summary. 



&/9 


a ''{nip) 


a- 'if.) 


./t / rrip 


fK/f. 




5.85 


1.49(29) 


1.62(27) 


0.154(40) 


1.206(15) 


14 


6.0 


1.89(14) 


2.06(15) 


0.155(11) 


1.205(15) 


10 


6.3 


2.96(25) 


3.24(27) 


0.155(19) 


1.14(14) 


10 



Using Eq. neglecting the anomalous contribution, and using the defini- 
tion of the decay constant, we can determine the pseudoscalar decay constant 
from the measurement of (0| jP). The results are summarized in Table 2. 
As with Bk, the estimates of the physical ratios fir/mp and Jk/ fir indi- 
cate good scaling behavior. They are also consistent with experiment, within 
rather large statistical errors. The errors in Table 2 are crude estimates de- 
rived using the statistical uncertainties only; the small sample size precludes 
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Fig. 2. The kaon B parameter. 

us from properly accounting for the correlations in the data. All of the results 
are so called effective values, calculated from the two-point correlators with- 
out fits and averaged over suitable plateaus. In the last column at the two 
larger couplings, an attempt has been made to carry out a more sophisticated 
jackknife error analysis all the way through to the final ratio. 

While all of the above results indicate good scaling, they must be checked 
further with improved statistics and a fully covariant fitting procedure. Also 
systematic effects like finite volume still need to be investigated. Only then 
can the continuum limit can be reliably taken. We note that a recent precise 
calculation using quenched Wilson quarks by the CP-PACS collaboration 
yields values of /,r and /k in the continuum limit that are inconsistent with 
experiment ||lO|| . 

In Fig. 3 we show the pion mass squared as a function of to. For TVs = 
10 the data at = 6.0 and 6.3 are consistent with chiral perturbation 

theory. However, at — 5.85, extrapolates to a positive non-zero 

value in the chiral limit for iV^ = 10 and 14. There is a large downward 
shift in the line as Ng goes from 10 to 14, but it still does not pass through 
the origin. However, at Ng = 18, to^ extrapolates to -0.004(19) at m = 0. 
The anomalous contribution on the r.h.s. of Eq. drops from 0.0098(5) to 
0.0038(2) as A''s varies from 10 to 18. It is interesting to note that at Ns = 10 
the anomalous piece is more than double the value at = 6.0 whereas the 
value at Ng = 18 is roughly the same. When the anomalous term is sufficiently 
small compared to the bare parameter m, the chiral symmetry is effectively 
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Fig. 3. The pion mass squared. Ns — 10 (octagons), 14 (squares), 18 (crosses). 

restored. From the above, it seems we must have {J^qln) / (JsItt) ^ 0.1 m. Of 
course, as m — > 0, one must also take Ns larger, which is analogous to the 
situation with the ordinary spatial volume. Fewer sites in the extra dimension 
may be sufficient at 5.85 if M is increased still further. 

Research supported by US DOE grant DE-AC0276CH0016. The numeri- 
cal computations were carried out on the NERSC T3E. 
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